      SUBROUTINE SETUP
      RETURN
      END
C   HS56
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS56'/
      DATA (X(I),I=1,7)/3*1.D0,3*.5097396789D0,0.9851107835D0/
      DATA N/7/ , NH/4/ , NG/0/
      DATA DEL0/1.00D0/ ,TAU0/1.D1/ ,TAU/.1D0/
      DATA ((GUNIT(I,J),I=1,3),J=0,4)/-1,0,0,-1,0,0,-1,0,0,
     F                                -1,0,0,-1,0,0/
      END
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX =-X(1)*X(2)*X(3)
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      INTEGER J
      DOUBLE PRECISION X(NX),GRADF(NX)
      ICGF=ICGF+1
      GRADF(1)=-X(2)*X(3)
      GRADF(2)=-X(1)*X(3)
      GRADF(3)=-X(1)*X(2)
      DO      100      J=4,7
      GRADF(J)=0.D0
  100 CONTINUE
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
       IF(I .EQ. 4)      GOTO 100
       HXI=X(I)-4.2D0*DSIN(X(I+3))**2
       RETURN
  100  CONTINUE
       HXI=X(1)+2.D0*X(2)+2.D0*X(3)-7.2D0*DSIN(X(7))**2
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      DO      10      J=1,7
      GRADHI(J)=0.D0
   10 CONTINUE
      IF(I .EQ. 4)      GOTO 100
      GRADHI(I)=1.D0
      GRADHI(I+3)=-8.4D0*DSIN(X(I+3))*DCOS(X(I+3))
      RETURN
  100 CONTINUE
      GRADHI(1)=1.D0
      GRADHI(2)=2.D0
      GRADHI(3)=2.D0
      GRADHI(7)=-14.4D0*DSIN(X(7))*DCOS(X(7))
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      RETURN
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I 
      DOUBLE PRECISION  X(NX) ,GRADGI(NX)
      RETURN
      END

